J 1^ production in nuclear collisions: measuring the transport coefficient 
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The observed pT-dependence of nuclear effects for J produced in heavy ion collisions at RHIC 
might look puzzling, since the nuclear suppression seems to fade at large pt- We explain this by 
the interplay of three mechanisms: (i) attenuation of Jj'^ in the hot medium created in the nuclear 
collision; (ii) initial state shadowing of charmed quarks and attenuation of a cc dipole propagating 
through the colliding nuclei; (iii) a strong Cronin effect for J/^ caused by saturation of gluons in 
the colliding nuclei. All three effects are well under control and calculated in a parameter free way, 
except for the transport coefficient go characterizing the medium. This is adjusted to the J/'^ data 
and found to be in good agreement with the pQCD prediction, but more than an order of magnitude 
smaller than what was extracted from jet quenching data within the energy loss scenario. 
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I. INTRODUCTION 

The recent measurements [H-O of J/^* produced in 
heavy ion collisions at RHIC have revealed unusual fea- 
tures of the transverse momentum distribution. While 
all species of hadrons measured so far demonstrate nu- 
clear suppression, which increases with pT and then levels 
off, the nuclei-to-pp ratio for J /"^ production, plotted in 
I rises with pT and has even a tendency to exceed 
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FIG. 1: (Color on line) Nuclear ratio Raa for central (0- 
20%) copper-copper (full circles and squares, upper curve) 
and gold-gold (empty circles, bottom curve) as function of 
transverse momentum of the J/^P. The curves are calculated 
with Eq. as is described in text. 



No explanation has been proposed so far, except for 
an exotic one [4] assuming that part of the production 
rate comes from accidental coalescence of c and c pairs 
available in the medium. Even if this might happen, one 
should consider first of all the conventional explanations, 
based on known dynamics. 

We consider here three different mechanisms affecting 
the production rate of J /"^ in heavy ion collisions: (i) fi- 
nal state attenuation (FSI) of J/^* in the dense medium; 
(ii) initial state interaction (ISI), nuclear shadowing of 
charm quarks and the breakup of the cc dipole propagat- 
ing through the colliding nuclei; (iii) ISI Cronin effect for 
J/^E" caused by gluon saturation in the colliding nuclei. 
All three effects certainly exist and are important, and 
below we present their evaluation, which is performed in 
a parameter free way, except the transport coefficient 
characterizing the hot medium. This is assumed to be 
unknown and is adjusted to reproduce the data with the 
value of go ~ 0.3 - O.SGeVVfm. This value IS an or- 
der or two of magnitude less than what was extracted so 
far from high-pT pion suppression observed in gold-gold 
collisions at ^/s = 200 GeV, and interpreted within the 
energy loss scenario [1, 0| • 

Only the first of the three effects mentioned above, the 
J/ 4' attenuation due to FSI, was considered in the recent 
publication [8], but the ISI suppression was ignored. Be- 
sides, the cc separation was assumed to be fixed during 
propagation through the medium, while the J /"^ wave 
function is fully formed within a very short distance, half 
a fermi (see next section III A| . As a result, go was grossly 
(5 times) overestimated in [8]. 

On the contrary, in [9] it was assumed that J /"^ is sup- 
pressed only by ISI, but propagates with no attenuation 
through the produced dense matter. The observed nu- 
clear effects were explained by ISI and by the suppressed 
feed-down from the decays of heavier states (x, 5'')j which 
can be dissolved in the hot medium. Such an approach 



does not look self-consistent: if J/'^ is absorbed even 
in the cold nuclear matter, it should be even more sup- 
pressed propagating through a dense medium. 



II. FINAL STATE INTERACTION OF J/* 
A. Time evolution of a small dipole 



A cc dipole is produced aX xp — in the cm. of 
the collisions, with a short time scale ~ \/ -\- p'^ 
and with a small transverse separation r ~ l/rric. Then it 
evolves its size and forms the J/ wave function. The full 
quantum-mechanical description of this process is based 
on the path integral technique |10|]. However, just a rough 
estimate of the formation time is sufficient here, since this 
time scale turns out to be very short. 

A small size dipole is expanding so fast that its initial 
size is quickly forgotten. Indeed, the speed of expansion 
of a dipole correlates with its size: the smaller the dipole 
is, the faster it is evolving. This is controlled by the 
uncertainty principle, k ~ l/r. 



B. J/^ attenuation in a dense medium 

A charmonium propagates a path length L in a 
medium with the survival probability 



S{L) = cxp 



- / dla[r{l)]p{l) 



(5) 



Here the path length and time are related as I = vt, 
with the J/^" speed v — \/l — (2mc/-E)^. The medium 
density is time dependent, and is assumed to dilute as 
p{t) = poto/t due to the longitudinal expansion. 

The dipole cross section for small dipoles is cr(r) = 
C , where r is the transverse cc separation. Corre- 
spondingly, CTj/^ = |C and we rely on the result 
of the realistic model [HI, [l^ for the mean J radius, 
\J {f'^j/x^) = 0.42 fm. The factor C for dipole-proton in- 
teractions is known from DIS data. Its value for a hot 
medium is unknown, as well as the medium properties. 
However, the factor C also controls broadening of a quark 
propagating through the medium [l^ [3] , 
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where E 
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- 2E* is the J/^* energy in the cm. of the 
collision; k ^ l/r is the transverse momentum of the c- 
quark relative to the J/^* direction. The solution of this 
equation reads 
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where Tq '-^ l/{p\-\-m1) is the initial dipole size squared, 
which is neglected in what follows. 

According to © the expanding cc reaches the J/vf size 
very fast. 



< 0.6 fm, 



(3) 



for J / 5* transverse momenta up to 5 GeV. This is about 
the expected time of creation of the medium. 

Another estimate of the formation time scale in the 
cm. of collision is llOll. 
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This 



where m^* is the mass of the first radial excitation, 
results in the same estimate as (jS)). 

We conclude that what is propagating through the 
medium is not a small cc dipole (pre-hadron), but a fully 
formed J/^*. 
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Thus, the factor C is related to the transport coefficient 
q [5i] , which is in- medium broadening per unit of length. 



C = 



2p 



(7) 



So one can represent the survival probability of J/^ in 
the medium, Eq. ([5]), as 



S{L) = exp 



3 V ,//*/ 



dlq{l) 



(8) 



The transport coefficient depends on the medium den- 
sity, which is a function of impact parameter and time. 
We rely on the conventional form 15|, 
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where b and r are the impact parameter of the collision 
and of the point where the q is defined. The transport co- 
efficient go corresponds to the maximum medium density 
produced at impact parameter t = in central collision 
(6 = 0) of two nuclei, at the time t — to after the colli- 
sion. In what follows we treat the transport coefficient % 
corresponding to the medium produced in central gold- 
gold collision at 6 = r = 0, as a adjusted parameter. It 
is rescaled for other nuclei according to the number of 
participants Upart (b, f) [l5l | . In what follows we consider 
collision of identical nuclei, A = B, at 6 = 0. 
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Eventually, integrating the attenuation factor Eq. ^ 
over different direction of propagation of the J/"^ pro- 
duced at impact parameter t one gets the FSI suppres- 
sion factor in the form, 



nFSI 



[t.Pt) = y -3(^7/*) j dlq{f+l) 

^ lo ^ 

(10) 

Here \ f+ = + P + 2tI cos 4>; and Iq = vto- The time 
scale to for creation and thermalization of the medium 
is rather uncertain, since gluons with different transverse 
momenta are radiated at different coherence times. We 
rely on the usual estimate to — 0.5 fm. 

The results are depicted by dotted curve in Fig. [5] for 
copper-copper collisions. We use qo = 0.45GeV^/fm, 



III. INITIAL STATE INTERACTION: 
SHADOWING AND ABSORPTION 

A. Higher twist shadowing of charm 

The same time scales, production and formation, look 
very different in the rest frame of one of the collision 
nuclei. While a fully formed J/^* propagates through 
the hot medium, in this reference frame a cc dipole with 
a size "frozen" by the Lorentz time dilation propagates 
through the cold nuclear matter. Both the production 
and formation times become longer by the Lorentz factor 
7 = 2i?,p/-\/4TOf"+"p^. The coherence time of cc pair 
production reads. 
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FIG. 2; Data for the nuclear ratio Raa for central copper- 
copper collisions as function of transverse momentum of the 
,7/^ from [2] (circles) and 0] (squares). The dotted curve 
shows the FSI effects, Eq. pop . The dashed curve includes the 
ISI effects, charm shadowing and absorption, Eq. p9p . The 
solid curve is also corrected for the Cronin effect Eq. ((28]). 

which allows to reproduce data well, provided that 
other corrections, discussed in the following sections, are 
added. 



and is rather long. At ^/s = 200 GeV it varies from 13 
to 4 fm for < < 5 GeV, i.e. is of the order of the 
nucleus size, or longer. Correspondingly, the formation 
time is even longer, ~ 5ic ^ Ra- 

A full calculation of the nuclear effects for J/'^ pro- 
duced in pA collisions, including the effects of shadowing 
and breakup interactions of the final cc, has not been 
done so far. Only production of the Xi the P-wave char- 
monium, which is a simpler case, was calculated in detail 
in [1^. Besides, a considerable fraction of J/^s are pro- 
duced via decay of heavier states, 5*', Xj etc. For our 
purposes it would be safer to use the experimental value 
of nuclear suppression observed in d — Au collisions at 
= 200 GeV 18]. Unfortunately the experimental un- 
certainty is still large, so we fix RdAi^F = 0) = 0.8 which 
is about the central value. 

Even if this nuclear suppression factor integrated over 
impact parameter is known, it is not sufficient to perform 
calculations for AA collisions. One has to know the b- 
dependence of Ra- Since no relevant data are available 
so far, we can only rely on the theory, being constrained 
by the integrated value of Ra- 

Since the coherence time Eq. (jlip in the rest frame of 
one of the colliding nuclei is rather long, we assume that 
the cc transverse separation is "frozen" by Lorentz time 
delation. This grossly simplifies the calculations. In the 
dipole approach the nuclear suppression factor caused by 
initial state c-quark shadowing and attenuation of the cc 
dipole, has the form [l(i |. 
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Here Ta{t,z) = J^^dz' pa{t, z'); e± = {er,±ey)/y/2 
is the polarization vector of the projectile gluon. The 
light-cone wave function of a cc in the gluon "^g^cdr) cc 
Ko{mcr) [l^, where Kq is the modified Bcssel function. 
Thus, the mean transverse size squared of a cc fluctuation 
of a gluon, (r^) = is small, at least an order of 

magnitude smaller than that of charmonia. Apparently, 
in a convolution of two r-distributions, narrow and wide, 
the mean size is controlled by the narrow one. 

The cross section adin(r) in (IT2t is the universal dipole- 
proton cross section [19|, which we use in the small-r 
approximation Udipij) « C{x2)r'^. The factor C{x2) is 
calculated in [23| as function of 0:2, which is the fractional 
light-cone momentum of the J/^P relative to the target. 

Shadowing of the process g ^ cc is controlled by the 
three-body g-q-q dipole cross section, which can be ex- 
pressed via the conventional dipole cross sections, 

9 1 

<^s{r) = - [(Jdip{ri) + (Jdip{r2)] - -^(Jdipiri ~ ^2), (13) 

where ri and r2 are the transverse vectors between the 
gluon and the q and q respectively. We neglect the dis- 
tribution of the fractional light-cone momentum of the q 
and q, fixing it at equal shares. Then ri = — r2 = so 
that 
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Here we also rely on the small-r approximation. 

Since we fixed the overall suppression RdA at the mea- 
sured value, and need to know only the impact parameter 
dependence, a rough estimate of Eq. ([T^ should be suf- 
ficient. Therefore, we approximate the result of integra- 
tion over r in ([T2|) replacing the dipole cross sections by 
an effective cross section, Udipir) =J> <yeff^ which we can 
adjust to reproduce RdA- Then the suppression factor 
Eq. (fT2)) takes the form. 
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The first term in square brackets represents shadowing, 
the second one is related to the survival probability of 
the produced colorless cc dipole. 

We extracted the value of (Jeff comparing the inte- 
grated suppression. 



Rna = J (frTAir) Rna{t), 



(16) 



with data [Hi for deuteron-gold collisions at = 
200 GeV, xp — 0. We fixed the measured ratio at 
RdAu = 0.8, and found (Jeff = 2.3 mb. 

This value can be compared with the theoretical ex- 
pectation. As was mentioned, in the convolution of the 
narrow distribution g^ccif) with the large size charmo- 
nium wave function, the latter can be fixed at r = 0, and 
the mean separation is fully controlled by the cc distri- 
bution in a gluon. Then the mean separation squared of 



a produced cc pair, i.e. a fluctuation which took part in 
the interaction, is given by 
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Now we are in a position to evaluate the effective cross 
section, 

(Jeff = C{E) {r'). (18) 

The energy dependent factor C{E) is calculated in [20| . 
At the energy of J/* E = 300 GeV (xp = 0, ^/s = 
200 GeV) this factor varies between C = 4.5 in the lead- 
ing order, down to C = 3.5 if higher order corrections 
are included. Correspondingly, the effective cross section 
Eq. ([T8|) range is 2.5 mb > Cef / > 2 mb, which is in excel- 
lent agreement with the value extracted from the RHIC 
data. 

We calculate the ISI suppression nucleus-nucleus colli- 
sions assuming that the suppression factors due to simul- 
taneous propagation of the cc pair through both nuclei 
factorize. We ignore the possible dynamics which can 
breakdown this assumption [21r,23j , so that the ISI sup- 
pression factor for a collision of nuclei A and B with 
impact parameter b reads. 



f (pTTA{T)TB{b- T)RNA{f)RNB{b- f) 
Jd^TTA{T)TB{b~T) 



, (19) 

Thus, the initial state interactions cause the additional 
suppression, Eq. (jl9p . of J/^* produced in heavy ion 
collision. The combined effect of ISI and FSI suppres- 
sion in copper-copper central collision is shown by the 
dashed curve in Fig. While it agrees with the data at 
Pt < 3 GeV, there is indication that data at higher pT 
are underestimated. 



B. Leading twist gluon shadowing 

Besides quark shadowing, which is a higher twist effect 
and scales as I /ml, the leading twist gluon shadowing, 
which depends on me logarithmically, may be important, 
depending on kinematics. In terms of the Fock state 
decomposition, gluon shadowing is related to higher Fock 
components in the projectile gluon, e.g. g — >■ qqg. Even 
this lowest state is heavier than just a qq and should 
have a shorter coherence time. In terms of Bjorken x 
this means that the onset of gluon shadowing is shifted 
towards low X2 compared to quark shadowing. Indeed, 
calculations [13] show that no gluon shadowing is possible 
above X2 ~ 10~^. Moreover, it was found in [25] that 
the coherence length which controls the onset of gluon 
shadowing, is scale independent, i.e. it the same for light 
and heavy quarks. This result of ^25^ can be understood 
via the energy denominator for the g — >■ qqg transition 
amplitude, 

1 
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where k and ag are the transverse and fractional light- 
cone momenta of the radiated gluon, respectively. The 
factor Ug, which is predominantly small, suppresses the 
mass term in ^U\\ . Besides, the mean transverse mo- 
mentum of gluons was found in jU to be rather large, 
^ (fc2) — 0.7 GeV. This is dictated by data on large 
mass diffraction, which is strongly suppressed compared 
to usual pQCD expectations. This phenomenon has been 
known in the Regge phenomenology as smallness of the 
triple Pomeron vertex. The large value of a/ (fc^), also 
supported by many other experimental evidences [26|, 
leads to suppression of the gluon radiation amplitude 
Eq. (PO)) and weak gluon shadowing [23|. The latter 
is confirmed by a NLO analysis of the DIS data [28| . 
but contradicts the recent analysis of [29^ , which resulted 
in a very strong gluon shadowing breaking the unitarity 
bound [30| . 

The coherence length available for gluon shadowing 
can also be related to the value of X2 [25[, 



which leads to a new more correct value X2 ~ 2x2. With 
the corrected minimal value X2i'>] = 2) = 0.0 05 g luon 
shadowing correction is tiny, just a few percent [24 133 1. 



IV. BROADENING OF GLUONS, CRONIN 
EFFECT 

In pA collisions projectile gluons propagating through 
the nucleus experience transverse momentum broaden- 
ing due to multiple collisions. As a result, the mean 
transverse momentum of produced charmonia is larger 
than in pp collisions. The dipole approach [ij, [13] is 
rather successful predicting broadening for heavy quarko- 
nia [13] and heavy Drell-Yan dileptons [131, = 
{Pt)pA — {Pt)pp ^ parameter free way, relying on the 
phenomenological cross section [13] fitted to photopro- 
duction and DIS data. Broadening for a gluon of energy 
E propagating a nuclear thickness Ta reads [20| . 



P9{l-Xi) 

X2 niM 
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Here the denominator presents the usual loffe time scale, 
which is as long as (0.2fm)/a;2, and may exceed the nu- 
clear size at small X2- The factor 1 — xi is usually ne- 
glected, assuming that xi is small, which is not always 
the case. More important is the factor « 0.1 eval- 
uated in [IE]- Its smallness is actually due to the large 
intrinsic transverse momentum of gluons, which we have 
just discussed above. 
Usually X2 is defined as 



X2 



V/* 



(22) 



It varies with pseudorapidity and reaches a minimum at 
the largest measured value of 77. At = 0, with the mean 
value of {p^) = 4GeV^ [ll], one gets X2 = 0.02, which is 
certainly too large for gluon shadowing. Therefore we can 
safely disregard this correction in further calculations, 
done at ?7 = 0. 

Notice that of course X2 decreases with 77 and reaches 
its minimal value X2 = 2.5 x 10~^ at the maximal ra- 
pidity ij = 2. Although this value of X2 allows some 
amount of gluon shadowing, we expect a tiny correction. 
Indeed, within the color singlet model (GSM) 3l] and its 
modified version [33 | , which provides so far the only suc- 
cessful description of J/^* production in pp collisions, the 
actual X2 is considerably larger than the value given by 
the usual definition Eq. . This is because in the GSM 
J/ 4* is produced accompanied by a gluon, and their total 
invariant mass Mgj/^ is considerably larger than mj/i^,. 



With the mass distribution, da/dM'^ 
gets the mean invariant mass squared. 
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where the proton saturation scale is 

/ E 



QqN{E) = 0.19 GcV : 
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0.14 



(24) 



(25) 



In fact, the broadening Eq. (|24|) is the saturation scale 
in the nucleus calculated in the leading order, i.e. without 
corrections for gluon saturation in the medium. Those 
corrections lead to about 20% reduction of (2G] . 

Remarkably, broadening does not alter the shape of 
the pT-distribution of produced J/^*. Indeed, data on 
pp, pA and even AA collisions, at the energies of fixed 
target experiments 35] and at RHIG [ll|, are described 
well by the simple parametrization. 
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where (jp^) is the mean transverse momentum squared, 
which varies dependent on the process. Therefore, the 
simplest way to calculate the pT-dependence of the nu- 
clear cross section would be just making a shift A in 
the mean value (p^) for pA compared to pp, where A is 
broadening given by Eq. (|24)) . Then the nuclear ratio as 
function of pt gets reads. 
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where RpA is the pA over pp ratio of pT-integrated cross 
sections, Eq. (|T6| . 

This simple procedure looks natural, although is not 
really proven. We can test it with the precise data from 
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the E866 experiment at Fermilab at ^/s = 39 GeV. All 
the input parameters in (j27l) are known from the same 
measurements (sgI \3T\ and other experiments at the same 
energy [38l - l40| and also from our calculations, Eq. 
(p^) = 1.5 GeV^; A = 0.08 GeV^ x A^/^; RpA = A-o o^. 
The A-dependence of the nuclear ratio calculated with 
Eq. ([77)) as function of px is compared with E866 data in 
Fig.[3]for the exponent characterizing the A-dependence, 
a — I + ln(i?pA)/ This comparison confirms the 
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FIG. 3: The exponent a — 1 + ln{RpA) / \n{A) as function of 
Pt calculated with Eq. (|27[1 in comparison with data from the 
E866 experiment [3^ . 

validity of the chosen procedure, at least within the mea- 
sure interval of pt ■ 

Inspired by the good agreement with pA data, we apply 
the same procedure to evaluation of the Cronin effect for 
J/^ production in central AA collisions. One can use for 
Raa{pt) the same Eq. (I?7l) . which should be modified 
replacing A and (p^) by corresponding values for A A 
collisions at ^/s — 200 GeV. Also one should replace the 
ratio of the pT-hrtegrated cross sections RpA R^A^i^ — 
0), which is calculated with Eq. ([T9l) . 



A = 0.93 GeV for gold-gold collisions. Notice that the 
mean value of transverse momentum squared in pp colli- 
sions at = 200 GeV, {p\) = 4GeV^, is considerably 
larger than at y/s = 39 GeV. 

Eventually, we are in a position to combine all three 
effects and calculate the pT-dependence of the nuclear 
ratio in central AA collisions, 



i?l//(5 = 0,pT) 



/ dVTl(r) R^'J{t,pt) RTI{r,PT) 
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(29) 

The result is depicted by solid curve in Fig. [2] in compar- 
ison with data for copper-copper collision. Calculations 
and data for gold-gold collisions are also shown in Fig. [T] 
Notice that the procedure Eq. ([77)) has not been con- 
fronted with data above pt = 4 GeV, so our extrapola- 
tion and predicted steep rise of the ratio at px, which 
might create problems with fcT-factorization, is not well 
justified. For this reason we tried another way to imple- 
ment broadening into the pT-distribution. An alternative 
way would be a simple convolution of broadening, which 
we take in the Gaussian form, with the px distribution 
in pp collisions. Then the nuclear modification factor 
R^AAi''''PT) get the form. 
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The results for copper-copper and lead-lead collisions are 
plotted in Fig. We see that the description of data at 
Pt < 5 GeV is unchanged compared to what was depicted 
in Fig. [TJ This means that our determination of go from 
the data is stable against the choice of the way how the 
Pt broadening is included. Only the behavior at larger 
Pt, where available data have poor accuracy, is altered, 
showing a weaker Cronin enhancement. 
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where the ratio of py-integrated cross sections, R^j^J{b — 
0, t), is given by Eq. (|19p without integration over r. 

According to Eqs. ([M|) - ([^5|) broadening slowly rises 
with energy. However, the J/^ energy in the nuclear 
rest frame is about the same in the E866 experiment 
{{E) = 230 GeV) and in Phenix data at xp = and 
Vs = 200 GeV {E = 330 GeV), so we neglect the differ- 
ence. Thus, in A A collisions A simply doubles compared 
to pA, and we get A = 0.64 GeV^ for coper-coper and 



V. PROBING DENSE MATTER AT SPS 

The nuclear suppression caused by FSI of the J/^ with 
the produced medium was determined in the NA50 and 
NA60 experiments comparing the measured nuclear sup- 
pression Raa with what one could expect as the cold nu- 
clear effects in initial state extrapolating from pA data. 
The latest results from the NA60 experiment [4l| for 
maximal number of participants corresponding to cen- 
tral collisions show that the FSI suppression factor is 
0.75±0.7. This experimental uncertainty is shown by the 
horizontal stripe in Fig. [Sj The curve shows dependence 
of the FSI suppression factor i?^^^(6 — 0), calculated 
with Eq. ([2^ . on the transport coefficient. The factor 



6 




2 4 6 8 

Pt (GeV) 



10 



FIG. 4: The same as in Fig. [T] but the Cronin effect is calcu- 
lated differently, with Eq. ((SO)! . 



jjISI 



was excluded from (I29p and integration over pT per- 
formed. Together with the experimental uncertainty this 
curve provides an interval of values of = 0.23 ± 0.07 
(vertical stripe), which allows to describe the observed 
"anomalous" suppression. It is about twice as small as 
we got from RHIC data. 
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FIG. 5: The FSI attenuation factor for J/^ produced in cen- 
tral lead-lead collisions at SPS. The curve corresponds to sup- 
pression versus (Jo calculated with Eq. (1291 with excluded ISI 
factor and integrated over pr- The horizontal stripe shows 
the magnitude and uncertainty of suppression reported by 
the NA60 experiment [4l[. The vertical stripe shows the in- 
terval of values of go which allow to describe the observed 
suppression. 



VI. SUMMARY AND OUTLOOK 

We performed an analysis of data for pT-dependent nu- 
clear effects in J/^* production in central copper-copper 
and gold- gold collisions observed at RHIC. These data, 
looking puzzling at first glance, have not received a 
proper interpretation so far. We evaluated the final state 
attenuation of the produced J/'^ in the created dense 
medium relating it to the transport coefficient, i.e. broad- 
ening of partons propagating through the medium. The 
key point, which allows to establish this relation, is the 
dipole description of broadening [l^, ■ 

The observed nuclear effects in J/^P production in AA 
collisions is interpreted as a combination of FSI of the 
fully formed J/^ in the dense medium, and the ISI ef- 
fects in production of J /"^ caused by multiple interac- 
tions of the colliding nuclei. The latter includes, attenu- 
ation of the produced cc dipole propagating through both 
nuclei, higher twist shadowing of charm quarks, and lead- 
ing twist gluon shadowing. Besides, gluon saturation in 
nuclei leads to a considerable broadening of gluons, which 
causes a strong Cronin effect for J/^*. This explains the 
observed unusual rise of the nuclear ratio with pT- 

All effects are evaluated in a parameter free way, except 
for the unknown properties of the produced hot medium. 
We employed the conventional model for the space-time 
development of the produced matter relating it to the 
number of participants. The only parameter adjusted to 
data, go, is the transport coefficient corresponding to a 
maximal density of the matter produced in central gold- 
gold collision. We found that the J/^E" data from RHIC 
are well reproduced with go ~ 0.3 — 0.5GeV^/fm. This 
is close to the expected value go — 0.5GeV^/fm and 
more than order of magnitude less than was found from 
jet quenching data within the energy loss scenario [Tt'I. 

We also examined the J/^* data from the NA60 experi- 
ments at SPS, which are available for py-integrated cross 
sections, and with already separated ISI effects. From 
the observed suppression in central lead-lead collisions 
we found go ~ 0.23 ± 0.07 GeVV fm. 

We performed the calculations assuming direct J/'^ 
production, but it is known that about 40% comes from 
the feed-down by decays of heavier charmonium states, 
X and VP' [i^l . Those states are about twice as big as the 
J/vP [l^l, and correspondingly should have a larger ab- 
sorption cross section. Therefore, adding these channels 
of J/ 4' production can only reduce the value of go, i.e. 
the above values should be treated as an upper bound. 
The bottom bound can be estimated assuming that the 
r^-approximation is valid up to the size of these exci- 
tations (which is certainly an exaggeration). Then the 
bottom bound for go is the above values times a factor 
0.7. 

We conclude that pr-dependent nuclear effects for J/^ 
production in heavy ion collisions provide a sensitive 
probe for the dense medium produced in these collisions. 
Both experimental and theoretical developments need 
further progress. More accurate pA (or dA) data are 
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needed for a better control of the ISI effects in AA col- 
lisions. Data for other heavy quarkonia T) would 
bring new precious information. 

On the theoretical side, more calculations are required 
to describe the observed centrality and rapidity depen- 
dence of nuclear effects. The azimuthal asymmetry can 
also be calculated. The small-r approximation for the 
dipole cross section used here may not be sufficiently ac- 
curate for a very dense matter. One should rely on a more 
elaborated r-dependence. A more realistic model for the 
space-time development of the dense medium, including 



transverse expansion, should be considered. 
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